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The existence of the light gluino - a color octet, Majorana fermion, of mass a few GeV 
or less, which is the superpartner of the gluon, is the subject of intensive discussion in the 
literature The possible impact of light gluinos on LEP/SLC experiments has been 

discussed recently @.[|] . 

One of the best ways to look for a trace of a still undiscovered light gluino is to evaluate 
its contributions to the precisely measured quantities, such as the hadronic decay widths of 
the Z bozon and the r lepton, using the perturbation theory. Perturbative QCD (pQCD) 
calculations of the above quantities up to 0(a^) have been completed in recent years 0- 
0. For a recent review of the status of these calculations see [[TO 
contributions see the review 



For the Electroweak 
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In the present paper we report results of the evaluation of pQCD corrections up to 0(a^) 
to the above quantities due to the gluino. We assume that the gluino is light and its mass is 
anywhere from less than 1 GeV up to several GeV and the squarks are so heavy that they 
decouple for the energy range considered (~ Mz). Note however that near the r mass, the 
gluino mass has to be well below M T , in order our results to be valid for the r decay rate. 

We start with a brief outline of the theoretical structure of the Z boson total hadronic 
decay rate to 0(a^) 
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Here the summation index runs over light quark flavors / = u,d, s, c, b. We define Xf = 
m 2 f(M z ) I 'M| and X t = m 2 (M z )/Mz, where we use the MS definition of quark masses. The 
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vector and axial couplings of quark / to the Z boson are Vf 
The electroweak self-energy and vertex corrections are absorbed in the factors pf and kf. 
The status report for the electroweak contributions is given in ref. jD]]. The Y^ A (Xf) are 
vector and axial parts of the well known Born approximation of V z (see, e.g., 0]). The 
terms S^(^(a,Xf) represent contributions due to diagrams involving at least one photon 



exchange between the quarks (see, e.g., pT|). The QCD contributions are represented by so 
called nonsinglet <5qcd() an d singlet £ V/A () terms. The singlet part is due to Feynman graphs 
with the electroweak currents in separate quark loops mediated by gluonic states. The other 
type graphs form the nonsinglet contribution. Ng in the above equation is the number 
of light gluinos that can appear virtually in some topological types of graphs, starting at 
0(a 2 ). The above QCD terms are calculated up to O(a^) within the standard model, with 



no gluinos 



(see, e.g., jlQl ). In the present work we calculate the QCD corrections 



up to the four-loop level that involve light gluino contributions. 

The calculational methods are very similar to that for standard QCD (for a detailed 
description see |[L2|| ). The hadronic decay rate of the Z boson in the tree level approximation 
for the electroweak sector can be evaluated as the imaginary part 

Y z = -^- £ lmU(m f ,m t ,s + iO) (2) 
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where the function II is defined through a correlation function of two flavor diagonal quark 
currents 

i J d 4 xe^(Tjl(x)jl(0)) = g^(m f ,m t ,Q 2 )-Q, l Q u U\m f ,m t ,Q 2 ). (3) 

Here, Q 2 is a large Euclidean momentum ~ — M§. The standard neutral weak current of a 
quark / coupled to the Z boson is jf = (GFM z /2y/2) 1 ^ 2 (vf'qj-'j ll qf + a/fZ/Tyys?/)- Because 
of this structure of the neutral weak current, the H-function may be decomposed into vector 
and axial parts n(m/, rrit, Q 2 ) = Tl v (mf,m t ,Q 2 ) + Tl A (rrif,mt,Q 2 ). For further calculational 
convenience we use the approximation 

n^K, m t , Q 2 ) = Il^(Q 2 , log m t ) + 0$) + 0(%) + £ (4) 

This is a legitimate expansion since our problem scale is Mz- Thus we work in the limit 
of zero light quark mass and infinitely heavy top quark mass. The last term in the above 
equation is the nonperturbative contribution, parametrized by semi-phenomenological quan- 
tities - the so called vacuum condensates. C^ A are their coefficient functions, that can also 
be evaluated perturbatively (see, e.g., |12| and references therein). In this work we ignore 



these contributions. The corrections due to light quark masses, especially for b quark, and 
finite top mass are not negligible. They are known up to the order considered (see, e.g., 
|lCfl and references therein). Note that Uq^ a (Q 2 , logm t ) (below we omit the subscript and 
superscripts for IT) still depends on m t . This is due to a logarithmic dependence of the 
axial singlet part on the top mass that is not suppressed by inverse powers of m t . Thus the 
decoupling of heavy particles is not manifest in the MS type prescriptions. The known mass 
dependent corrections can simply be added to our result. On the other hand, in the limit of 
vanishing light quark masses, the vector and axial parts of the nonsinglet contributions are 
identical. We calculate first the nonsinglet vector (axial) part and then we treat the singlet 
axial part. 

The further steps seem to be straightforward. We write the diagrammatic representa- 
tion for II and calculate relevant multiloop Feynman graphs analytically using dimensional 
regularization [|l3j and the modified minimal subtraction (MS) prescription fll]. After the 



renormalization of coupling we can get the final result using eq.(|2|). Unfortunately this 
straightforward scheme is realistic only up to three-loop level, since even the advanced com- 
puter program HEPLoops [Tj|| for evaluation of multiloop Feynman graphs can evaluate 



diagrams only up to three-loops. Fortunately, with the aid of the renormalization group 
and a unique feature of the MS prescription one can reduce the four-loop calculation to the 
evaluation of only one-, two-, and three- loop graphs. The detailed outline of the method 
and further references are given in ref. Jl2| . 



The running strong coupling obeys the renormalization group equation a s /3(a s ) = 
fj^dag/dfi 2 , where the QCD f3 function coefficients for MS type schemes are defined as 
follows (3{a s ) = -AjK/4tt) - fcfa/Aic) 2 + 0{a 3 s ), (3 = llC A /3 - 4TN/3 - 2C A N- g /3, 
fa. = 3ACI/3 - 20C A TN/3 - AC F TN - 16C^JV s /3. In the above coefficients we include 
light gluino contributions. The corresponding terms are obtained from the analysis of all 



diagrams contributing to one- and two-loop /3-function coefficients and known results (16 
are confirmed. 
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(a) (b) (c) 



FIG. 1. Four-loop Feynman graphs with light gluino contributing in the nonsinglet (a,b) and 
singlet (c) parts of Tz- The dashed lines correspond to light gluino propagators. The solid and 
wave lines correspond to quark and gluon propagators. 



The evaluation of the three-loop approximation to the decay rate including the light 
gluino effect is fairly trivial, since there are only two three-loop graphs where the gluino can 
appear. We use the three-loop graph-by-graph results of and have confirmed three-loop 
numerical results given in ^j. Note that there are no gluino contributions at one or two- 
loop levels. At the four-loop level, the situation is more complex. Although we use still 
unpublished four-loop diagrammatic results from the first work in ||, it was necessary to 
reanalyze all 27 four-loop graphs where the light gluino can appear. These graphs contribute 
with different color weights in the case of gluino. There are several graphs that required 



recalculation with the help of the HEPLoops program [15|. Two of those graphs are shown in 
Fig. la and Fig. lb. In the previous standard QCD calculations ||, contributions from those 
and several other graphs were taken into account with the help of the full two-loop gluon 
propagator inserted into a two-loop graph. In the case of a light gluino it was necessary to 
treat these graphs individually, because of different color and symmetry weights. 

For the nonsinglet part of the four-loop QCD correction to the Z decay rate, including 
the light gluino contribution, we obtained the following MS analytical result 
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In the above expression the logarithmic contributions are summed up into the running con- 
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stant by taking /i 2 = Mf . Those contributions can trivially be restored using the renormal- 
ization group (see @]). Inserting the standard SU C (3) eigenvalues of the Casimir operators 
for the fundamental and adjoint representations Cp = 4/3, Ca = 3 and also T=l/2, we 
obtain the following numerical result 
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For the singlet part the light gluino contribution to 0(a^) shows up in only one single 
graph (Fig.lc) in the axial channel. In the vector channel, the corresponding graph vanishes 
due to Furry's theorem | 1~7| . The nonvanishing of the graphs like in Fig.lc is due to a large 
mass splitting in the top-bottom doublet. This effect was first evaluated exactly in the 
three-loop level || and then extended to four-loop using the large mass expansion method 
f|. We use the result of f| to extract the light gluino contribution to the four- loop axial 
singlet part. We obtain 
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Note that throughout this paper the strong coupling is defined for five flavor effective theory. 

For completeness, we also give the result for singlet vector part in the limit of vanishing 
light quark masses and infinitely heavy top mass , although the gluino does not contribute 
here 
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The terms of 0{X^ V ) || are about two orders of magnitude less than the leading term || 
and are completely negligible. 

Next, we use our result to calculate the light gluino contribution to the r-lepton decay 
rate to 0{a^ s ) in perturbative QCD. We consider the familiar ratio R T = r(r~ — » v T + 
hadrons)/r(r~ — > v T e~V e ). Here we are interested only in perturbation theory contributions 
and we do not consider nonperturbative and instanton corrections. For the calculational 
method and references see ||12|| . As before, we work in the limit of vanishing light quark 
masses and infinitely large heavy quark masses. Note that for the scale ~ M T , u,d,s quarks 
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are considered as light and c,b,t are heavy quarks. We use our diagrammatic results obtained 
for the Z boson case to evaluate correlator of the charged weak currents of quarks coupled 
to W boson. We obtain the following result for QCD perturbative contributions to R T , 
including a light gluino 
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and numerically we get 

R T (Ml) = 3(0.998 ± 0.002) 
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The three- and four-loop corrections due to a light gluino are very large and they increase 
the extracted a s (M T ) by more than 15%. (On the extraction of a s (M T ) from the r-lepton 
decay width see, e.g., fl8| and references therein.) For the case Ng — Eqs. (||) and (|i~0"D 
agree with the known results 0. 

Summarizing, we have calculated the light gluino contribution to hadronic decay rates of 
the Z boson and the r-lepton to four-loop level in perturbative QCD. The corrections in the 
case of the Z boson decrease the three- and four-loop coefficients by about 25% each. The 
net effect of a light gluino is to increase a s (M z ) by about 2%. In the case of the r-lepton, 
the corrections are very large and they decrease the three- and four-loop coefficients by 22% 
and 74% respectively. The r hadronic decay rate remains a major counter-indication to the 
hypothesis of a light gluino. The value of a s (M T ) extracted from r decay with a light gluino 
and extrapolated to Mz overestimates the direct measurement of a s (Mz)- Therefore, if a 
light gluino exists there must be appreciable, positive contributions to the r hadronic width 
from, for instance, the nonperturbative region or other as-yet-unknown sources. 
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Note Added After this paper was completed, we received a paper [HI], where the light 
gluino contribution was calculated to four-loops for the vector part of the Z boson decay 
rate. The result is in complete agreement with our eq.(H). 
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